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Sensing the internal dynamics of individual nuclear spins or clusters of nuclear spins has recently become 
possible by observing the coherence decay of a nearby electronic spin: the weak magnetic noise is amplified 
by a periodic, multi-pulse decoupling sequence. However, it remains challenging to robustly infer underlying 
atomic-scale structure from decoherence traces in all but the simplest cases. We introduce Floquet spectroscopy 
as a versatile paradigm for analysis of these experiments, and argue it offers a number of general advantages. In 
particular, this technique generalises to more complex situations, offering physical insight in regimes of many- 
body dynamics, strong coupling and pulses of finite duration. As there is no requirement for resonant driving, the 
proposed spectroscopic approach permits physical interpretation of striking, but overlooked, coherence decay 
features in terms of the form of the avoided crossings of the underlying quasienergy eigenspectrum. This is 
exemplified by a set of “diamond” shaped features arising for transverse-field scans in the case of single-spin 
sensing by NV-centers in diamond. We investigate also applications for donors in silicon showing that the 
resulting tunable interaction strengths offer highly promising future sensors. 

PACS numbers: 


I. INTRODUCTION 

There is enormous interest in the rapidly advancing field of 
detection and imaging at the single spin level JUS), mainly 
with NV colour centers but also other defects in diamond 
ElEI not only as a source of versatile qubits for quantum in¬ 
formation m m and entanglement generation fill , but prin¬ 
cipally because they underpin a new generation of quantum 
sensors, for magnetometry and atomic scale characterisation 
of the environment 112-18]. In the widely-studied case of dy¬ 
namical decoupling quantum sensing, a sequence of pulses 
modulates the coherent evolution of the sensor and in some 
cases, sharp “dips” in coherence allow one to detect, and infer 
useful characteristics of, nearby single spins or small spin- 
clusters. Complexities in the environment being studied mean 
that the single isolated sharp dip is found in a restricted sub¬ 
set of the data and motivates the development of more general 
or alternative methods of analysis. In particular, many decou¬ 
pling sequences are temporally periodic. 

In this case, Lloquet’s theorem provides a canonical form 
for the solution of periodically driven systems and has found 
wide applicability in various branches of quantum physics 
since 1965 CD, especially in light-matter interactions with 
continuous driving and multi-photon atomic physics. Llo¬ 
quet’s theorem can be applied to any periodic quantum Hamil¬ 
tonian for which H{t + r to t ) = H(t). In practical imple¬ 
mentations, instead of analysing the eigenstates of the static 
Hamiltonian, which are appropriate only in the perturbative 
limit of weak driving, one employs instead the eigenspectrum 
of the one period time-evolution propagator. Lloquet theory 
is employed in analysis also of Nuclear Magnetic Resonance 
(NMR) spectra and related (20lL where applications are essen¬ 
tially limited to resonant and continuous driving. 

But this approach has not been considered for analysis of 
coherence behaviors in this new generation of multipulse spin 
sensing experiments. We argue that the Lloquet spectroscopy 
method proposed is better adapted to regimes of strong quan¬ 


tum entanglement between the sensor and detected spin sys¬ 
tems than signal processing methods applied to classical ac 
signals; or geometric approaches based on two-state systems. 
In this work, we find that Lloquet theory can augment current 
methods of analysis in several ways: 

(1) Lloquet theory is equally applicable to off-resonant as 
to resonant driving. Understanding of spin-sensing data is of¬ 
ten cast in signal processing terms: the multi-pulse sequence 
imposes a filter function which selects an ac signal with a rea¬ 
sonably well-defined characteristic frequency cc ac /27r which 
may, in turn, be used to infer interatomic coupling parameters. 
When c o ac is resonant with the pulse interval r (see Lig.[ljc)), 
so when cc ac /7r = l/(4r), a narrow dip in coherence is ob¬ 
served. Here we show that away from such resonances, or 
even for broadened resonances, Lloquet theory can shed in¬ 
sight on other striking features which are not narrow dips 
but can, nevertheless, yield rich information for atomic scale 
characterisation. The key reason is that we show the widths 
and shapes of these features may be understood in terms of 
avoided crossings of an underlying Lloquet spectrum. 

(2) Current experiments operate primarily in regimes of 

weak-coupling. The pulses involve consecutive switching be¬ 
tween two electronic states u,d\ the associated characteristic 
frequencies of the detected spins cu n , ujd can vary significantly, 
since, for stronger coupling, there is significant back-action 
and entanglement between the sensor and detected spins. Lor 
weak coupling, uj u ~ c Od and in addition average Hamiltonian 
theory models apply, predicting typically uj ac « + Ud). 

Lloquet theory remains valid regardless of coupling strength; 
we examine regimes of failure, obtain alternative forms for 
uj ac and show that the avoided crossings shed insight in these 
regimes. It remains also valid even if there is non-trivial evo¬ 
lution due to finite duration of the pulses, a problem only re¬ 
cently identified fl7l . 

(3) Lor detection of two-state systems (spins or spin pairs 
which reduce to an effective pseudospin) geometric methods 
fT6l 2Tj] are widely used to interpret data and yield analytical 
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FIG. 1: (a) Current experiments have employed the S = 1 electronic spins of NV-centers to successfuly detect (i) single nuclear spins 

115: .21 |22) (ii) the internal dynamics of nuclear spin pairs m as well to characterise on the atomic scale, by estimating parameters such 
as electronic-nuclear dipolar couplings A and inter-nuclear dipolar couplings C 12 . (b) Generic sensor detecting a pair of nuclear spins: the 
electronic spin state is in a superposition of “up” | u) and “down” | d) states. The nuclear dynamics and its characteristic frequencies uo u ,d 
depend on the electronic state. In turn, the electronic coherence is sensitive to the resulting weak ac noise from the nuclei. This may be 
amplified by dynamical decoupling control such as CPMG, leading to observed “dips” in coherence. These are at well-defined frequencies in 
typical weak coupling regimes where the nuclear dynamics is not too different in the u , d subspaces. However, strong coupling regimes do not 
necessarily yield sharp dips, (c) Additional complexities occur for pulses of finite durations (d) It is is also challenging to differentiate between 
(i) independent pairs of spins and (ii) many-body effects from an equivalent interacting cluster. Floquet theory is not restricted to single spins 
or spin pairs and can be applied also to analysis of larger, correlated spin clusters, strong coupling and off-resonant driving. 


expressions for the coherence decays. Although here the Flo¬ 
quet method already sheds additional physical insight, its full 
value is that it is universally valid even for higher dimensional 
state-spaces so would facilitate studies of e.g. multi-spin clus¬ 
ters. 

The key-features of dynamical-decoupling based quantum 
sensing, using a multipulse periodic sequence, are illustrated 
in Fig|T] A 7t/2 pulse prepares the sensor system in a super¬ 
position state ^(t = 0) = ^(|?z) + |d)) 0 S(0), where B( 0) 
is the detected spin-cluster at initial time. In turn, interaction 
with the sensor means that the spin-cluster becomes entangled 
with the sensor ~ ^(\u)B(u u , t) + \d)B(uod, t))> where 

B(co u ,d, t) = (T U: d) Np B( 0), for a pulse sequence (with prop¬ 
agator T u ^) which is repeated N p times. The detected spin 
dynamics is associated with a characteristic frequency which 
is state-dependent. The temporal coherence C(t) = (5 + ), is 
given by C (t) = (B (uj u ,t)\B (ujd ,t) ) to within a normalisation 
factor; averaged over bath states, it simulates the experimen¬ 
tally measured signal. 

II. FLOQUET THEORY 

Floquet’s theorem is generally applicable to periodically- 
forced systems, classical or quantum, but it allows one specif¬ 
ically to write solutions to the Schrodinger equation in terms 
quasi-energy states (QES) \ipi(t)) = exp (— ieit)\$i) where 
e t is the quasi energy, | $;(£)) = | $i(t + r tot )), r tot is the 
period and l = 1,Z2 (Z3> is the dimension). However, for 
problems (such as our present study) where we require only 
“stroboscopic” knowledge of our system (i.e. read-out once 
every period T tot ) , the solution is even simpler. We can obtain 
Floquet phases/modes simply as the eigenvalues/eigenstates 


of the one-period unitary evolution operator T(r t01 ,0). The 
Floquet modes obey the eigenvalue equation: 

T(r tot )|^) «A Z |$ Z > =exp(—i£ z )|$ z ) (1) 

where now Si is the eigenphase (the Floquet phase) and e/ = 
Si/r tot is the quasienergy. For sensing, we can obtain Flo¬ 
quet phases/modes simply as the eigenvalues/eigenstates of 
T U: d, the basic periodic sequence; for example, for the CPMG 
sequence in Fighl r t01 = 4r. In that instance Ti\$u) = 

e~' l£ i | $>u) where i = u,d denotes the state of the sensor 
spin. The eigenphases for N p = 1 fully determine the time- 
evolution of the system since if the pulse sequence is repeated 
N p times, we just scale the eigenphases, so for longer time- 
propagation: 

(£)*” ( 2 ) 

In the present work, we show that these Floquet eigen¬ 
phases and eigenstates have particular important properties: 

(a) The eigenvalues are the same for the upper and lower 
states, i.e. e lS i = e lS i = e l£l . This holds even for pulses 
of finite duration, in typical cases. In other words, the evolu¬ 
tion of the Tu Np ^B( 0) and T^ Np ^B{ 0) are characterised by the 
same set of effective frequencies e/ = Si/r tot , in contrast to 
typical static, geometric approaches where two distinct sets of 
frequencies uj u i and ujdi are involved. 

(b) The eigenvectors do not, in general, coincide but we show 
that (to within a phase term), the eigenvectors are related to 
each other by a half-period evolution e.g. T u (r to t/^)\^ui} 

(c) Minima in coherence (of prime importance for sensing, 
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FIG. 2: (a) usual Geometric approach: under CPMG-N control ((b)) the detected spins represent two-state systems which precess about 

effective magnetic field, depending on the “up” | u) and “down” | d) states of the probe spin. The coherence dips are understood by following 
the precessions and relative angles between these spins, with increasing N. (c) Spectroscopic picture. The dips in coherence occur at avoided 
crossings of the Floquet eigenstates. Both the position and contrast of the decoherence dip is related to the curvature of the crossing. This is 
characterised by the splitting between the states 25 and the deviation from the early r 0 linear evolution. The early time evolution (the eo 
quasienergy) gives the the dip position t<h v for average Hamiltonian theory; the coherence minimum is given by £(^ R ) — (d) The dip 

contrast depends on the degree of curvature of the crossing, characterised by the level-repulsion strength, 5 — 2 _ £(^p)^. (e) 

NY-center decoherence “diamonds”. While typical experimental studies scan along parallel field (u z ) component (thus remaining in weak- 
coupling, single-dip regime), scanning the transverse magnetic field ( uj x ) would produce diamond pattern of high decoherence regions, as 
avoided crossings widen (and even overlap) then narrow (here uj z = 0 and A\\ = 50kHz). Upper panel shows full oscillating coherence 
function, for N p — 10 pulse pairs, lower panel shows coherence envelopes, filled as N p —)• oo. Here uj z = 0. Boundaries of the diamonds 
trace out (green) r — 7r/2 (ujd + w u ) and (cyan) r = 7r/2 (ujd — cj u ) (see below), (f)) expanded version of low field region showing shape of 
avoided crossings versus coherence traces corresponding to two cuts (i),(ii), indicated in the upper panel. 


whether sharp dips or not) occur at avoided crossings of Flo¬ 
quet eigenstates, where e l£l ~ e lSk . Once the Floquet phases 
and modes are obtained, one can obtain the general form of 
the decoherence for arbitrary times which, averaged over bath 
states, yields: 

1 D 

(C(t = N p 4r)> = pJ2 e ~ iNpi£l ~ £l,) \Zdi'\$ui) I 2 (3) 
1 , 1 ' 

Derivations of (a)-(c) are given in the Appendix. Although 
properties (a)-(c) are quite generic, physical insight on the 
Floquet picture is more easily gained from two-state sys¬ 
tems, where direct comparison with usual geometric methods 
[2JJ, 23| 24] is also possible. For the two-state case, eigen¬ 
values must be conjugate pairs A± = e ±l£ . Level cross¬ 
ings occur when A + ~ A_ hence the crossings must occur 
at £ ~ 0,7T, 2tt. The generic properties of states at avoided 
crossings (see Appendix) then imply coherence dips occur at 
£ ~ 7 r. We now first investigate the Floquet dynamics for 
these two-state single-spin or single pseudospin models. 


model, noting that for single-spin detection there is a gen¬ 
uine spin, while for pair-dynamics |14j 16] it is a pseu¬ 
dospin. It has led to a successful, widely used geometric 
model (see Fig|2ja) and (34l ) where the evolution of the pseu¬ 
dospin is conditional on the state i = it, d of the probe and 
corresponds to precession about an effective magnetic field: 
Hi = |lT • a = \[Xg x + Zi<j z ) where cr X: cr z are Pauli ma¬ 
trices in the usual spin basis; in the pseudospin case of course, 
we have | t|) I t) and | It) I !))• The X , Zi de¬ 
pend on the physical system (see lf34l for details); but for 
NV centers h n ~ (cu x ,0, Ay + u z ) while h d ~ (ju x , 0, uo z ) 
where = (uj x , 0,u; z ) is the external magnetic field, and 
A|| the parallel component of the hyperfine interaction. For 
spin pair-sensing, on the other hand, IT = ^(C 12 ,0, AAi) 
where AAi = 2(Ai — A2)(i\S z \i) represents the energy de¬ 
tuning between the nuclear spins in the pair and S represents 
the operator for the sensor spin. The eigenvalues of Hi are 

l o U:d = ±^^Jx 2 + d and the orientation of the effective 
field is 6i = arctan(X/z^). For two state systems, we obtain: 


A. Single spin or spin pair detection 

Both pair flip-flop dynamics as well as single spin- 
dynamics (in systems like NV centres where a crystal field 
leads to non trivial one-spin dynamics) can be approximated 
by a two-state Hamiltonian. We term this the pseudospin 


(At)) = 1-2 


'cos 2 [£(t)/2] - cos 2 [£ (t/2)] " 
cos 2 [£{ t )/ 2 ] 


sin 2 [N p £(t)] 
(4) 


This is a key result of the present work as it means one can 
give the full coherence function using only the Floquet phases. 
To calculate the Floquet eigenphase £(r), as well as its half- 
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FIG. 3: Coherence decay behavior for an electron spin detecting a 
flip-flopping pair of nuclear spins, for a donor in silicon system (see 
(Ml) with tunable interactions. C(Bo,t) exhibits a rich structure 
in the two-dimensional r, Bo plane which is not evident in the nor¬ 
mal traces at constant Bo. Decoherence map is shown for different 
R = AA/C12 , (A A = Ai — A 2 ). Large R corresponds to weaker 
dipolar coupling C 12 and the maps trace the locus of a set of isolated 
sharp dips in coherence. For smaller R, there are no longer single 
dips; nevertheless the envelopes (given by F(r)) are well-defined 
and track the behavior of the underlying Floquet avoided crossings. 
The background oscillatory structure depends on N p , the envelopes 
do not. Time t = 4 N p r\ (colour scale linear, with black = 1, yel¬ 
low < 0.5). Similar behavior is obtained for several transitions of 
Si:Bi and other donors, but specific parameters are for 12 —)► 9 ESR 
transition of Si:Bi and 2N P m 40. 


period value £{r / 2) in Eqp] one may use cos(£ (s)) = 
cos(2c d u s) cos(2c dds) ~ sin(2a; u 8) sin(2cu^ < s) cos (0 U — 0 d ) 
with s = r ors = r/2. Thus the coherence takes the form 
C(t) = 1 — 2 F(t) sin 2 [N p £ (r)], which is the product of a 
smooth envelope F(r), independent of N p , superimposed on 
a fast oscillating function sin 2 [N p £(r)], dependent on N p . 

Full comparison with geometric methods are in lf34l where 
we argue it is the condition 


£(r dip /2) = tt/2 (5) 

which best specifies the dip positions. The depth of the dip is 
related to the eigenvalue splitting parameter S = n — £(r d i p )', 
at the dip 

C(t = 4 NpT dip ) = 1-2 sin 2 (N p 5). (6) 

Hence, if £ = 7 r, there is no dip, so a true level crossing 
provides no signal. For N p 5 > tt/2 the width and shape of 
the dip becomes largely independent of N p and is fully de¬ 
termined by the Floquet anti-crossing and envelope function, 



Time = 4 N p r, (ms) Time = 4 N p r, (ms) 


FIG. 4: Comparisons between the dip positions obtained with Eq@ 
(blue line) and average Hamiltonian theory Eq[7] (red line) for the 
full coherence function (left panels) as well as its envelope (right). 
Within the field sweep there are global weak-coupling points (eg 
Bo ~ 0.19 T) where there is weak coupling regardless of the clus¬ 
ter properties and where the decoherence envelope collapses into a 
single sharp dip, a useful feature if high resolution is required: here 
there is always good agreement with average Hamiltonian theory. 
These points correspond to so called optimal working points |30 » 31 j 
of silicon donors. Hence, the advantage of such systems as future 
spin sensors, in addition to their very long ~ 1 s coherence times, is 
that a magnetic field sweep could tune the dynamics from the weak 
to strong coupling regimes, (a) and (d) R=100 (b) and (e) R=20 (c) 
and (f) R=10. 


since the sin 2 (N P S) prefactor simply superposes fast oscilla¬ 
tions on F(t). A narrow avoided-crossing (low splitting, S 
small) gives a single, sharp (but weaker) coherence dip, while 
a large S crossing has a broad envelope. It is only for low 
N P S <C 7t/2 that the dip height is strongly dependent on N p ; 
here the central height increases as ( N P S ) 2 . 

Comparison with Average Hamiltonian models (see t34l 
for details). A frequently used approximation in spin sens¬ 
ing is the average Hamiltonian model whereby the eigenval¬ 
ues cuav of, \ {H U + Hd), the time-averaged Hamiltonian pro¬ 
vide an estimate of the dip positions and that T d i p = 

From Fig|2](c) and (d) we can equate linear behaviour in our 
eigenvalues (narrow crossing, linear shape) both with the oc¬ 
currence of a sharp dip as well as validity of the averaged 
Hamiltonian model. In particular, for small r, £{r) 4eor 
corresponds to the averaged Hamiltonian results. Expanding 
the cos£(t) from below Eq|4], for small r, we easily obtain 
e 0 = \(w 2 +u 2 d + 2 Lo u u d cos (i 9 U - Od)) 1 / 2 and thus: 

Tdtp 2(uj 2 + uj + 2 uj u u d cos (9 U - Od)) 1 / 2 ^ 

Expressing quantities in terms of the pseudofield components 
X, Z u , d we can show that this is equivalent to the expression 
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FIG. 5: Fingerprinting multiple environmental spin cluster-pairs via their decoherence“bar-codes” illustrates the effect of 3-body correlations. 
The figure shows the coherence as a function of magnetic field B 0 and pulse interval r, calculated with a full numerical propagation under the 
total Hamiltonian for N p = 100. (a) denotes three independent pairs while (b)shows three interacting spins, with otherwise equivalent dipolar 
couplings and intrabath interactions as illustrated in Fig[^b). One evident difference (and signature of a cluster of three spins) are the doublets 
due to the two separate subspaces of the three interacting spins. The splittings are directly related to the interactions. For the 3-cluster, in fact 
there is a secular contribution from interactions between spins, greatly amplifying their contribution. The two right hand panels (i) and (ii) 
show single traces corresponding to the cuts in (b) as well as the six corresponding eigenphases: in case (i) in a weak coupling regime the 
dips are narrow and the eigenphases behave like three independent pairs; the eigenvalues correspond to conjugate pairs (with blue, red and 
green lines denoting the three pairs). In case (ii) there is stronger coupling, the avoided crossings of the corresponding eigenphases are broader 
giving raise to the features shown in the coherence maps. 


cc av = ^((X 2 + ( z u+ z d ) 2 ^i /2 use( j j n S pi n _detection experi¬ 
ments [14, 16] and to cc av — \{w u + u d ) f° r 6u — Od¬ 
in. APPLICATIONS 

A. Experimental control of quasienergy crossings 

The above motivates us to investigate possibilities for ex¬ 
perimental control of the avoided crossings, by varying S, 
even in the simple one-spin or spin-pair case. In typical sens¬ 
ing with NV-centers, we have cu z > uj x A, thus in Eq|7] 
we have 0 U ~ 0 d <C i r, thus r dip ~ 2 (aJ+c^) ■ However, 
setting uj z = 0 and increasing uo x causes the anti- crossings to 
widen and narrow successively, forming a checkerboard pat¬ 
tern of diamonds. This behavior was illustrated in Fig|2] (e) 
and (f). In particular Fig|2](e) illustrates the usefulness of the 
2D map; it is not easy to clearly discern the behavior from an 


individual trace (as in Fig [2] (f)) especially if N p is not very 
large. We note that higher harmonics have larger 5 than lower 
harmonics at the same parameters. 

However, here we consider in addition *5 = 1/2 systems as 
potential sensors. These might include silicon vacancies but 
in particular we focus on electron donors in silicon. Although 
techniques analogous to optical read-out and polarization of 
NV centres are not fully developed, there has been consider¬ 
able progress in single-spin detection f 251127 1 . These systems 
benefit also from extremely long coherence times (of order 
seconds) for cryogenically cooled samples. They are also an 
ideal test-bed for the theory as one can vary 0 Ul 0 d over a wide 
range as magnetic field B 0 is swept. For donor systems, the 
surrounding 29 Si nuclear spin dynamics does not generate an 
ac signal as there is no internal crystal axis, in contrast to the 
case of NV centers in diamond where the surrounding nuclear 
spins precess around an effective quantisation axis which is 
no longer only the external magnetic field. For donors, sin- 
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gle strongly coupled 29 Si nuclei have recently been detected 
via the static shift of the donor frequency 128]. However, the 
interesting coherence dynamics in these systems involve only 
pairs or larger clusters of spins where the flip-flopping dynam¬ 
ics generates an ac signal [29, 30]. 

Formally, the state-conditional dynamics for donors is very 
similar to that for NV centers: The dynamics correspond to 
an effective spin precessing about effective magnetic fields 
h z = !(Ci 2 , 0 , AAi) where A Ai = (A\ — A 2 )Pi ; But in 
contrast to NV centers Pi(Bo) = 2(i\S z \i), the polarisation 
of the state (see 1341 ) varies strongly with the magnetic field 
I3T1I321L while for NVs, (i\S z \i) = 0 , ±1 is fixed for the mod¬ 
est fields used in experiments. 

Fig [3] shows the field-dependence of the coherence for a va¬ 
riety of coupling strengths. The behavior may be compared 
with the NV-centers: in this case, the coherence dips trace out 
a curved locus, like NVs, for stronger X = C\ 2 /2 component 
of the pseudofield, the envelopes broaden, but there is a sim¬ 
ilar pattern of intermittent broadening and narrowing. There 
is a striking feature at (Bq = 188mT in Fig[3]where all deco¬ 
herence envelopes “collapse” to a sharp dip. This is one of a 
set of special fields (Optimal working points) where 0 U ~ 0 d 
and uj u ~ uod and which have been investigated theoretically 
and experimentally for their favourable coherence properties 

GolCm ED. 

But, in the present work, we find also that these points 
correspond also to very narrow Floquet avoided crossings, at 
which 5 —0. Fig |4] compares dip position predicted by aver¬ 
age Hamiltonian theory Eq|7] with the accurate dip condition 
£{r dip /2) = tt/2 . In Fig[4j A A = (Ai — A 2 ), (which approx¬ 
imately sets the timescale for R 1 ) was fixed, while C± 2 
(the intra-bath dipolar coupling) was varied to obtain different 
values of R = AAjC\ 2 . 

By means of a detailed theoretical analysis we can show 
that average Hamiltonian theory is valid if (i)|h n + h d \ ;$> 
|h n — h d \ or if (ii) \h u - h d | > |h n + h d |. Condition (i) cor¬ 
responds to the weak-coupling regimes typical of NV sensing 
experiments, where uo u ~ c o d and 0 U ~ 6 U \ it is also the 
regime of the Optimal Working Points, where average Hamil¬ 
tonian theory is always valid. Regime (ii) is not typical of 
sensing experiments; for the spin-1/2 system of Fig. [5] it 
would correspond to the spins nearly antialigned, thus P u ~ 
—P d . For Fig. ^condition (ii) implies (A A) 2 (P U — P d ) 2 > 
(C12) 2 + (A Ay(P u + P d ) 2 . In particular, for P u « -P d 
we obtain the condition | A A(P U — P d ) | \C± 2 |. Noting that 
| P u — P d | — 0.2 — 2 , this means that average Hamiltonian the¬ 
ory is valid at all fields for large R > 100 as seen in Fig. |4ja). 
However, for increasing intra-bath dipolar coupling C\ 2 , the 
theory ceases to be valid away from the small OWP regime in 
the center, as seen in Figs. |4](b) and (c) for smaller R. 

B. Detection of multi-spin clusters 

In this section we apply the Floquet approach to the system 
depicted in FigjTJd): we compare the decoherence “finger¬ 
print” of three independent spin pairs (analogous, formally, 
to the detection of three independent spins by NMR) with a 3- 


cluster which, in the absence of many-body interactions would 
give a similar signature. 

For the 3-cluster , we take three spins, with hyperfine cou¬ 
plings Ak = Ai,A 2: As to the sensor spin and with mutual 
dipolar interactions Cij = C12, C23, C31. Disregarding inter¬ 
actions, the energy cost of the spin flips is A^ = Ai — Aj. 



0 50 100 150 200 

time (ms) 


FIG. 6: Decoherence for an interacting cluster of three spins (3- 
cluster) . The coloured lines show comparisons with Eq[l0| showing 
excellent agreement with numerics obtained by diagonalisation of 
the full joint sensor-cluster Hamiltonian. 

For the independent pairs , we take three spin pairs, with 
the same dipolar interactions as the 3-cluster, but which 
are independent of each other. To have similar frequencies as 
the 3-cluster, we must have similar energy cost of a all three 
spin flips; and they must obey the cyclic condition of the 3- 
cluster A 12 + A 23 + A 31 = 0. Pair 1 has two spins with 
interaction C\ 2 and a pair of hyperfine couplings (Ai,A 2 ); 
pair 2 has interaction C 23 and hyperfine couplings (4 2 , A$)\ 
pair 3 has C 31 and hyperfine couplings (A 3 , Ai). We take 
C 12 = C 23 = kHz and C 31 = kHz, realistic 
values for nuclear impurities in the silicon lattice. We take 
A 1 = ^ kHz, As = kHz and A 2 = 0, thus our 
pairs correspond to R ~ 100 — 40 (as defined in Fig[3]) so 
the interactions are sufficiently weak to make their detection 
challenging but sufficiently strong to, below, illustrate impor¬ 
tant features. The choice of A 2 = 0 does not involve much 
loss of generality. If a state-dependent Hamiltonian is chosen, 
the A \, A 2 , As values can be shifted by an arbitrary constant 
without perturbing the dynamics. If the full Hamiltonian is 
considered, there can be higher order effects such as hyper¬ 
fine mediated corrections to C ZJ . This correction is very small 
for our parameters but is tested by full numerics below. 

Solution of total Hamiltonian First we set aside all pseu¬ 
dospin approximations and do diagonalisations of the full 
Hamiltonian followed by time-propagation, using the com¬ 
plete 8 -state basis of the 3-cluster as well as the complete 
basis of the bismuth sensor including the host nuclear spin. 
Thus, unlike Fig [3] we do not use the analytical form for the 
parameter P, \ it emerges from the numerics. We evaluate the 
decoherence numerically rather than using Eq|4] A similar 
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calculation was carried out with the three disjoint pairs, then 
the decoherence was averaged over the thermal ensemble of 
nuclear spin states (of which there are eight for the 3-cluster). 
Figure [5] shows maps of the coherence in the (r, Bq )-space in 
both cases. 

One conclusion to be drawn from comparisons between full 
numerics and the analytical (one-pair) Eq|4] is that the struc¬ 
ture in Fig[3] is surprisingly robust; without bath state aver¬ 
aging, full numerics give similar structure to Fig [3] (obtained 
from Eq|4]for one bath state). 

One striking feature of the 3-cluster decoherence map in 
Fig|5jb) is that some lines are split into “doublets” with very 
similar structure. The origin of these is in the average over the 
bath states; examining maps for the individual 8 bath states, 
we see that while the I z = ±3/2 cluster states | ttt) and 
| HI) make no appreciable contribution, the doublets arise 
from the separate I z = ±1/2 subspaces, which do not mix. 
In other words, the | ttl)> I tit) and | |tt) states with total 
quantum number I z = ±1/2 interact only weakly with the 
equivalent I z = —1/2 subspace, but each provides a locus 
of dips with a slightly different shift. In contrast to the spin 
pairs, in the case of the 3 -cluster, the secular Ising (Cjkhj Izk) 
components yield a non-trivial dynamical effect. 

Fig|5ji) and (ii) also shows a cut for two field values and 
compares with the behavior of the Floquet eigenstates. We can 
see that near the “weak-coupling” regime of optimal working 
points (i), the dips are sharp and narrow as are the avoided 
crossings; in contrast, away from the OWP point, avoided 
crossings are broader and even overlap (ii); the level splitting 
is much larger. We can estimate the point where two eigen¬ 
values will collide, and hence Tdi p by exploiting the fact that 
in either case, the average Hamiltonian theory value is not too 
far from the accurate value Tdip — Tdip • We estimate Floquet 
quasi-energies by considering only the diagonals. We obtain: 

e i = ^ — Aj — A k )(P u ± Pd) ± Cij ± Cik — Cj k ( 8 ) 


where z,j, k = 1,2,3 or cyclic permutations give ei= 1^,3 
quasienergies. 

Thus we estimate the dip positions from the fact that the 
quasi-energies represent the gradients of the spectral lines in 
Fig|5](i) and (ii), hence estimate the degeneracy point: 


(Im) 
' dip 


27r 


e l ~ e m 


(9) 


for the dip arising from the difference between the l and ra-th 
quasienergy. One finding is that the secular contribution from 
the dipolar coupling greatly amplifies the effect of the (usually 
weaker) Cij dipolar coupling between the nuclei, as it is a 
linear contribution. This is in contrast to disjoint pairs; if the 
dipolar coupling is weak, since oji — ± \ \J C\ 2 ± (Pi A A) 2 , 
for C 12 < Pi A A the non-secular contributions in the disjoint 
pairs represent a very small quadratic shift. 

In terms of the interaction strengths, the two dips of the first 
doublet correspond to: 

r ( ±12 ) ~_?!!_ do) 

dip |A 12 (P„ + P d ) ± 2(C 3 i — C 23 )| 


and similarly for other doublets. 

In Fig [6] we compare values from Eq[K)]with the full numer¬ 
ics. Thus the mean position exposes the value of A i2 while 
the splittings expose the dipolar coefficients. 


IV. CONCLUSIONS 


The extension of technologies such as MRI and NMR to 
the nano-scale is an outstanding technical challenge which is 
leading not only rapid experimental progress, but also the de¬ 
velopment of new methods to analyse data and to optimise 
information gathering on the atomic scale structure. 

Motivated by this, in the present work we introduce Floquet 
spectroscopy as an insightful new paradigm for understand¬ 
ing and analysis of spin sensing experiments. The approach 
is universally valid for any type of periodic driving whether 
resonant or not. Hence, here potential applications have been 
explored for analysis of different physical regimes and sen¬ 
sors which are not necessarily associated with a single sharp 
resonant ‘dip’ but may nevertheless potentially still offer well 
delineated features. 

Our key findings are (i) that there is an underlying struc¬ 
ture associated with Floquet avoided crossings and the Flo¬ 
quet spectrum which is potentially information rich; it repre¬ 
sents an envelope on the usually studied coherence dips with 
a shape controlled by the widths of the avoided crossings, 
(ii) that the Floquet approach clarifies also regimes where 
the commonly-used average Hamiltonian theory methods will 
fail, (iii) The method’s generality extends beyond single spin 
and pseudospin systems and is also useful for higher dimen¬ 
sional spin systems, and potentially any dynamical decoupling 
protocol, provided it is temporally periodic. 

Acknowledgements We are very grateful to Setrak Balian, 
Gary Wolfowicz and Gavin Morley for helpful discussions. 


Appendix A: Floquet spectrum 


1. Coherence minima and avoided crossings 

A key result of the present work is that the coherence 
dips associated with single spin sensing are associated with 
avoided crossings of the underlying Floquet spectrum and in 
this appendix, this conjecture is justified. 

In our study, we consider the important class of spin¬ 
sensing experiments for which an electronic sensor spin S 
is coupled to each environmental nucleus via the effective 
Hamiltonian: 


H ~ (i\S z \i)A • I (Al) 

Where i = u,d and A is a vector representing the hyperfine 
interaction. The dependence on S z only arises because of the 
large energy difference between electronic and nuclear states; 
in the case of NV-centers, the above is valid only for magnetic 
fields of magnitude and orientation which does not mix the 
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electronic states. The result is a state conditional Hamiltonian: 

H = -^=[| u)(u\ (8) H u + \d)(d\ (8) H d ] (A2) 
v 2 

where the Hi are the effective bath Hamiltonians discussed 
in Section [HA] in the main text so that an initial joint sensor- 
target spin state -0(0) ~ ^(\u) + \d))B(0) evolves into an 
(in general) entangled state: 

ip(t) = 2= (| u )B u (t) + | d)B d (t)). (A3) 

Experiments probe the coherence C(t) = (£ + ). While ex¬ 
perimental comparison involves averaging over thermal bath 
states Tr[pS + ], without loss of generality we consider a pure 
state £ = (ip(t)\S + \ip(t)) oc (B u (t)\B d (t)). 

Maximum entanglement occurs whenever (B u (t)\B d (t)) = 
0. However, the general condition for a minimum or dip to be 
observed is in fact \B u (t)) = — | B d (t)). This is regardless of 
the particular dynamical decoupling sequence applied. The 
key question for design of an experimental pulse sequence 
(say, CPMG) is for which pulse interval r and total pulse num¬ 
ber number N = 2 N p will correspond to underlying quantum 
evolution 

(B u (t = m p r)\B d (t = 4 N p r)} = -1 (A4) 

and thus a minimum in the function £, which to within an 
unimportant normalisation factor we take, C = (B u (t)\B d (t)) 
(we note that for decoherence experiments probing \C\, this in 
fact corresponds to a maximum in the coherence). 

The Floquet approach is based on the premise that for any 
periodically driven quantum system, the Floquet states 4>j ful¬ 
fil the same role as eigenstates of a Hamiltonian in a time- 
independent system. Thus if the initial quantum state is pro¬ 
jected into a Floquet basis, i.e. 6(0) = then its 

temporal evolution is known for all time. 

A central finding for the present work is that for the pulsed 
dynamical decoupling, the eigenspectrum is independent of 
the sensor spin state thus £^ = = £i(r) where the Flo¬ 

quet eigenspectrum £i(r) is a function of the experimentally 
chosen pulse interval r. 

The above results are proved in the next section, but we 
use them now to explain why coherence minima are associ¬ 
ated with avoided crossings. Since the Floquet spectra are the 
same, if there is an avoided crossing and thus a near degener¬ 
acy, e £l ~ e £k in the u subspace of sensor states, there 
will simultaneously be an avoided crossing in the lower d sub¬ 
space of sensor states. 

2. Avoided crossings for two-level system 

Although the eigenspectra are the same, in general the cor¬ 
responding eigenstates or Floquet states are not. 4>^(r) ^ 
<h u (r) for arbitrary r. Hence the temporal evolution: 

B u {t) = a u+ ^ u+ e~ iN ^ + a u -$ u -e+ iN ” £ ^ ± 
B d (t ) = a d+ ^ d+ e~ iN ^ + a d ^ d .e +iN ^ (A5) 


and thus entanglement with the sensor is established since the 
sensor-target B u , d (t) state is no longer factorisable. 

One exception occurs for r = 0, where all the Floquet 
states reduce to the unperturbed (thermal states). For a two- 
state system, without loss of generality, 4> u+ (r = 0) = 
4>d + (r = 0) = | t) or alternatively 4> w _(r = 0) = $> d -(r = 

0)U). 

Another, most interesting, exception is at a level crossing, 
where the eigenstates take the same form. It is a textbook 
result for level crossings (also known as Fandau-Zener transi¬ 
tions) that the the unperturbed states are maximally mixed and 
become sums and differences of the unperturbed states. The 
implication for the present case, is that the Floquet states for 
both upper and lower state must coincide at approximately 
^=(| \) ± | t))- This allows for two distinct possibilities: 

(i) In the first case, 

= $ d+ ~ -k(| |) + 11» and 

= *d- ^ T(| _ | t )) (A6) 

This possibility is the trivial case where the Floquet states 
for upper and lower sensor state are identical. There is never 
any difference in the evolution, no entanglement and so 
(B u (t)\B d (t)) = +1. 

(ii) In the second case, 

$„+ = ~ 4=(| + | t» and 

a>„_ = $ d+ ~-l (U) _| t » (A7) 

In this case, (B u \B d ) = cos2 N p £(r) which may for an ap¬ 
propriate choice of iVp «tt/ 2£ attain the minimal value for a 
dip (B u \B d ) = -1. 

Diagonalisation of the two dimensional unitary matrix is 
straightforward (see l34l ) and it is clear its eigenvalues must 
be conjugate pairs A± = e ±lS . For a two-level case, the 
avoided crossing condition is A + = A_ and hence coherence 
dips must occur at £ ~ 0,7r, 2tt with case (i) occurring at 
£ ~ 7r and case (ii) occurring at £ ~ 0, 2tt. 

3. Coherence minima for the general case 

For a spin cluster of arbitrary size, an initial pure state can 
be projected into the upper or lower Floquet basis: 

D 

|8(0)> = £<*"|S(0)>|<&«> (A8) 

l 

where |T>^) are Floquet eigenstates for upper {i = u) and 
lower (i = d) sensor state respectively and D is the dimension 
of the one-period unitary evolution operator. The state after a 
time t = 4 N p r is then 

D 

I =J2e- iN ^(^u\B(0mu) (A9) 

l 
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where Si are the eigenphases. 

Thus the coherence decay of the sensor spin (C(t) = 
(Bi(t)\B u (t))) is given in a Floquet basis by: 


C(t = N p 4r) 


D 

E* 

1 , 1 ' 


—iN p (£i —£ d ) 


h| 2 |($ d H^)l 2 , 


(A 10) 


Here \ai\ 2 = |($ Mi |i3(0))| 2 . 

The target spins are in fact in a thermal ensemble, which 
given small nuclear energy scale are all equally likely. Thus 
we must average over the thermal bath states and calcu¬ 
late (C) = (1 /D)Y^j £j> where Cj is the coherence 

function evaluated for the bath intitially in the j-th ther¬ 
mal bath state, \Bj(0)). Under the thermal average we find 

= 1. This produces Eq. [3] in the 
main text. While C for a pure state is complex, the (in any case 
small) imaginary part vanishes under the bath average and we 
consider only the real part of C. 

Returning briefly to simplest case of D = 2, treated in the 
previous subsection, Eq. ( |A10 ) may be rewritten as 

£i-£ 2 n 


(C(T)) = l-2\(*d-\*u+)f*n 2 K 


(All) 


which is of the form C(r) = 1 — F(r)f(N p ,r) where 
/(A/p, r) G [0,1] is a pulse number dependent oscillation. If 
we disregard the oscillation, we obtain the pulse number in¬ 
dependent minimal bound of the coherence function given by 
C em (r) = 1 — F(r), which we call the coherence envelope. 

Even if these envelopes are not necessarily sharp “dips” (es¬ 
pecially in strong-coupling regimes) they can correspond to 
sharply delineated structures (for both NV centers and donors) 
which should still be observable experimentally and can yield 
valuable information about the atomic-scale structure. We 
have shown above that the dips occur when |4> w+ ) — | 

For a general bath, of arbitrary dimension D , we can re¬ 
arrange Eq.( |A10| ) using only orthonormality and complete¬ 
ness of the eigenstates J2 k , \(&dk' 10)1 2 = 1 into paired con¬ 
tributions. 


( £ > = i - E[i( $<B 'i $ui >r+ k $ «hI $ «*')i 2 ] 

l<l' 

x sin 2 (N p — -F- ) (A12) 

This is, again, composed of pulse number independent 
envelopes superimposed with pulse number dependent 
oscillations. For a minimum, we require the term in square 
brackets to be maximised; this will occur at a level crossing 
between a given pair of eigenstates l and l' as argued in the 
previous subsection. For D > 2, level crossings between 
Si m Sir occur at arbitrary Si and no longer at Sip ^ it. 
Dips occurring at the point for which \(&di'\&ui}\ 2 = 1 and 
\(&di\&ui '}\ 2 = 1, generalises, to arbitrary dimension, the 
two-state orthogonality condition that \& u +) = |3>d-)- 

Appendix B: Symmetry of eigenphases for CPMG control 

For decoupling sequences like CPMG, the Floquet phases 
are independent of the sensor spin state, regardless of the di¬ 
mensionality of the bath states i.e. S = s[ d ^ = Si. To 


show this we first construct the basic propagator, for total pe¬ 
riod T to t = 4r, which is to be repeated periodically: 

f ( (2) (4r) = T {u) (r)f (d) (r)f (d) (r)f (tl) (r) = f {u) (2r)f {d) (2r) 

f ( (2) (4r) = T {d) (t)T { u) (t)T { u) (r)f (d) (r) = f (d) (2r)f (u) (2r) 

(Bl) 


~ (2) 

We can then obtain the eigenvalues for : 

f ( (2 ) ) (4r)|$ Mi ) = f (u) (2r)f {d) (2r)|$ uI > = e~ i£ ‘\$ ul ). 

(B2) 

Here exp (—iSi) is the Ith eigenvalue of If we apply the 
half period operator, T( d ) (2r): 


T(d) (2t )T( w) (2r )T(d } (2r) | <& u i) = e- <f *f (d) (2r)|$ ttZ ), 

(B3) 

this is equivalent to 

T ( (2) (4r)T (d) (2r)|^) = e~^f (d) { 2r)|$ uJ ), (B4) 


implies that T^(2r)\^ u i) is a eigenstate of , i.e. 

T^d){2r)\^ui) oc | $di}> where the factor of proportionality 
is a complex phase exp (ipid)- Similarly, 

T( u )(2r)|$dz) = ex.p(ipi u )\$ u i) (B5) 

for which pid + p,i u = Si. From this we see that each Floquet 
state | $ui), | &di) is the half-period evolution of the other (up 
to a complex phase). 


Thus exp (—iSi) is also an eigenvalue of T^. 


Eq 


B4 


Appendix C: Pulses of finite duration 

Provided the (4 r) can be decomposed into products 

of sub-propagators, as in Eq |C2| the it pulses do not have to 
be of very short duration. If we write: 

T (n) (2r + 25) = f iu) (r)T n (25)f {d) (r) 

T {d) (2r + 25) = f {d) (r)T n (25)f iu) (r) (Cl) 

Then we can write the full propagator in the same form as 
previously: 


T ( (2 }(4t') = T (u) (2r')T (d) (2rO 

= T(d) (2r , )T( u) (2t) (C2) 

but with r' = r + 5. Then all the above follow; the pulses 
of finite duration are still assumed to be it pulses but there 
can be some arbitrary evolution of the system during the finite 
interval 5, but properties such as the state independence of the 
Floquet phases still hold. 
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